Effect of Cholesterol on Biomimetic Membrane Curvature and Coronavirus Fusion Peptide Encapsulation.

Biomimetic cubic phases can be used for protein encapsulation in a variety of applications such as biosensors and drug delivery. Cubic phases with a high concentration of cholesterol and phospholipids were obtained herein. It is shown that the cubic phase structure can be maintained with a higher concentration of biomimetic membrane additives than has been reported previously. Opposing effects on the curvature of the membrane were observed upon the addition of phospholipids and cholesterol. Furthermore, the coronavirus fusion peptide significantly increased the negative curvature of the biomimetic membrane with cholesterol. We show that the viral fusion peptide can undergo structural changes leading to the formation of hydrophobic α-helices that insert into the lipid bilayer. This is of high importance, as a fusion peptide that induces increased negative curvature as shown by the formation of inverse hexagonal phases allows for greater contact area between two membranes, which is required for viral fusion to occur. The cytotoxicity assay showed that the toxicity toward HeLa cells was dramatically decreased when the cholesterol or peptide level in the nanoparticles increased. This suggests that the addition of cholesterol can improve the biocompatibility of the cubic phase nanoparticles, making them safer for use in biomedical applications. As the results, this work improves the potential for the biomedical end-use applications of the nonlamellar lipid nanoparticles and shows the need of systematic formulation studies due to the complex interplay of all components.

Comprehensive classification of proteins based on structures that engage lipids by COMPOSEL.

Structures can now be predicted for any protein using programs like AlphaFold and Rosetta, which rely on a foundation of experimentally determined structures of architecturally diverse proteins. The accuracy of such artificial intelligence and machine learning (AI/ML) approaches benefits from the specification of restraints which assist in navigating the universe of folds to converge on models most representative of a given protein's physiological structure. This is especially pertinent for membrane proteins, with structures and functions that depend on their presence in lipid bilayers. Structures of proteins in their membrane environments could conceivably be predicted from AI/ML approaches with user-specificized parameters that describe each element of the architecture of a membrane protein accompanied by its lipid environment. We propose the Classification Of Membrane Proteins based On Structures Engaging Lipids (COMPOSEL), which builds on existing nomenclature types for monotopic, bitopic, polytopic and peripheral membrane proteins as well as lipids. Functional and regulatory elements are also defined in the scripts, as shown with membrane fusing synaptotagmins, multidomain PDZD8 and Protrudin proteins that recognize phosphoinositide (PI) lipids, the intrinsically disordered MARCKS protein, caveolins, the ß barrel assembly machine (BAM), an adhesion G-protein coupled receptor (aGPCR) and two lipid modifying enzymes - diacylglycerol kinase DGKε and fatty aldehyde dehydrogenase FALDH. This demonstrates how COMPOSEL communicates lipid interactivity as well as signaling mechanisms and binding of metabolites, drug molecules, polypeptides or nucleic acids to describe the operations of any protein. Moreover COMPOSEL can be scaled to express how genomes encode membrane structures and how our organs are infiltrated by pathogens such as SARS-CoV-2.

Ionizable lipids penetrate phospholipid bilayers with high phase transition temperatures: perspectives from free energy calculations.

The efficacies of modern gene-therapies strongly depend on their contents. At the same time the most potent formulations might not contain the best compounds. In this work we investigated the effect of phospholipids and their saturation on the binding ability of (6Z,9Z,28Z,31Z)-heptatriacont-6,9,28,31-tetraene-19-yl 4-(dimethylamino) butanoate (DLin-MC3-DMA) to model membranes at the neutral pH. We discovered that DLin-MC3-DMA has affinity to the most saturated monocomponent lipid bilayer 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and an aversion to the unsaturated one 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The preference to a certain membrane was also well-correlated to the phase transition temperatures of phospholipid bilayers, and to their structural and dynamical properties. Additionally, in the case of the presence of DLin-MC3-DMA in the membrane with DOPC the ionizable lipid penetrated it, which indicates possible synergistic effects. Comparisons with other ionizable lipids were performed using a model lipid bilayer of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC). Particularly, the lipids heptadecan-9-yl 8-[2-hydroxyethyl-(6-oxo-6-undecoxyhexyl)amino]octanoate (SM-102) and [(4-hydroxybutyl) azanediyl] di(hexane-6,1-diyl) bis(2-hexyldecanoate) (ALC-0315) from modern mRNA-vaccines against COVID-19 were investigated and force fields parameters were derived for those new lipids. It was discovered that ALC-0315 binds strongest to the membrane, while DLin-MC3-DMA is not able to reside in the bilayer center. The ability to penetrate the membrane POPC by SM-102 and ALC-0315 can be related to their saturation, comparing to DLin-MC3-DMA.

The Cytoplasmic Domain of the SARS-CoV-2 Envelope Protein Assembles into a ß-Sheet Bundle in Lipid Bilayers.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) envelope (E) protein forms a pentameric ion channel in the lipid membrane of the endoplasmic reticulum Golgi intermediate compartment (ERGIC) of the infected cell. The cytoplasmic domain of E interacts with host proteins to cause virus pathogenicity and may also mediate virus assembly and budding. To understand the structural basis of these functions, here we investigate the conformation and dynamics of an E protein construct (residues 8-65) that encompasses the transmembrane domain and the majority of the cytoplasmic domain using solid-state NMR. 13C and 15N chemical shifts indicate that the cytoplasmic domain adopts a ß-sheet-rich conformation that contains three ß-strands separated by turns. The five subunits associate into an umbrella-shaped bundle that is attached to the transmembrane helices by a disordered loop. Water-edited NMR spectra indicate that the third ß-strand at the C terminus of the protein is well hydrated, indicating that it is at the surface of the ß-bundle. The structure of the cytoplasmic domain cannot be uniquely determined from the inter-residue correlations obtained here due to ambiguities in distinguishing intermolecular and intramolecular contacts for a compact pentameric assembly of this small domain. Instead, we present four structural topologies that are consistent with the measured inter-residue contacts. These data indicate that the cytoplasmic domain of the SARS-CoV-2 E protein has a strong propensity to adopt ß-sheet conformations when the protein is present at high concentrations in lipid bilayers. The equilibrium between the ß-strand conformation and the previously reported α-helical conformation may underlie the multiple functions of E in the host cell and in the virion.

Oligomerization-Dependent Beta-Structure Formation in SARS-CoV-2 Envelope Protein.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the current COVID-19 pandemic. In SARS-CoV-2, the channel-forming envelope (E) protein is almost identical to the E protein in SARS-CoV, and both share an identical α-helical channel-forming domain. Structures for the latter are available in both detergent and lipid membranes. However, models of the extramembrane domains have only been obtained from solution NMR in detergents, and show no ß-strands, in contrast to secondary-structure predictions. Herein, we have studied the conformation of purified SARS-CoV-2 E protein in lipid bilayers that mimic the composition of ER-Golgi intermediate compartment (ERGIC) membranes. The full-length E protein at high protein-to-lipid ratios produced a clear shoulder at 1635 cm-1, consistent with the ß-structure, but this was absent when the E protein was diluted, which instead showed a band at around 1688 cm-1, usually assigned to ß-turns. The results were similar with a mixture of POPC:POPG (2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine/3-glycerol) and also when using an E-truncated form (residues 8-65). However, the latter only showed ß-structure formation at the highest concentration tested, while having a weaker oligomerization tendency in detergents than in full-length E protein. Therefore, we conclude that E monomer-monomer interaction triggers formation of the ß-structure from an undefined structure (possibly ß-turns) in at least about 15 residues located at the C-terminal extramembrane domain. Due to its proximity to the channel, this ß-structure domain could modulate channel activity or modify membrane structure at the time of virion formation inside the cell.

SARS-CoV-2 Envelope Protein Forms Clustered Pentamers in Lipid Bilayers.

The SARS-CoV-2 envelope (E) protein is a viroporin associated with the acute respiratory symptoms of COVID-19. E forms cation-selective ion channels that assemble in the lipid membrane of the endoplasmic reticulum Golgi intermediate compartment. The channel activity of E is linked to the inflammatory response of the host cell to the virus. Like many viroporins, E is thought to oligomerize with a well-defined stoichiometry. However, attempts to determine the E stoichiometry have led to inconclusive results and suggested mixtures of oligomers whose exact nature might vary with the detergent used. Here, we employ 19F solid-state nuclear magnetic resonance and the centerband-only detection of exchange (CODEX) technique to determine the oligomeric number of E's transmembrane domain (ETM) in lipid bilayers. The CODEX equilibrium value, which corresponds to the inverse of the oligomeric number, indicates that ETM assembles into pentamers in lipid bilayers, without any detectable fraction of low-molecular-weight oligomers. Unexpectedly, at high peptide concentrations and in the presence of the lipid phosphatidylinositol, the CODEX data indicate that more than five 19F spins are within a detectable distance of about 2 nm, suggesting that the ETM pentamers cluster in the lipid bilayer. Monte Carlo simulations that take into account peptide-peptide and peptide-lipid interactions yielded pentamer clusters that reproduced the CODEX data. This supramolecular organization is likely important for E-mediated virus assembly and budding and for the channel function of the protein.

Continuum Models of Membrane Fusion: Evolution of the Theory.

Starting from fertilization, through tissue growth, hormone secretion, synaptic transmission, and sometimes morbid events of carcinogenesis and viral infections, membrane fusion regulates the whole life of high organisms. Despite that, a lot of fusion processes still lack well-established models and even a list of main actors. A merger of membranes requires their topological rearrangements controlled by elastic properties of a lipid bilayer. That is why continuum models based on theories of membrane elasticity are actively applied for the construction of physical models of membrane fusion. Started from the view on the membrane as a structureless film with postulated geometry of fusion intermediates, they developed along with experimental and computational techniques to a powerful tool for prediction of the whole process with molecular accuracy. In the present review, focusing on fusion processes occurring in eukaryotic cells, we scrutinize the history of these models, their evolution and complication, as well as open questions and remaining theoretical problems. We show that modern approaches in this field allow continuum models of membrane fusion to stand shoulder to shoulder with molecular dynamics simulations, and provide the deepest understanding of this process in multiple biological systems.

Strikingly Different Roles of SARS-CoV-2 Fusion Peptides Uncovered by Neutron Scattering.

Coronavirus disease-2019 (COVID-19), a potentially lethal respiratory illness caused by the coronavirus SARS-CoV-2, emerged in the end of 2019 and has since spread aggressively across the globe. A thorough understanding of the molecular mechanisms of cellular infection by coronaviruses is therefore of utmost importance. A critical stage in infection is the fusion between viral and host membranes. Here, we present a detailed investigation of the role of selected SARS-CoV-2 Spike fusion peptides, and the influence of calcium and cholesterol, in this fusion process. Structural information from specular neutron reflectometry and small angle neutron scattering, complemented by dynamics information from quasi-elastic and spin-echo neutron spectroscopy, revealed strikingly different functions encoded in the Spike fusion domain. Calcium drives the N-terminal of the Spike fusion domain to fully cross the host plasma membrane. Removing calcium, however, reorients the peptide back to the lipid leaflet closest to the virus, leading to significant changes in lipid fluidity and rigidity. In conjunction with other regions of the fusion domain, which are also positioned to bridge and dehydrate viral and host membranes, the molecular events leading to cell entry by SARS-CoV-2 are proposed.

Role of Ionizable Lipids in SARS-CoV-2 Vaccines As Revealed by Molecular Dynamics Simulations: From Membrane Structure to Interaction with mRNA Fragments.

Recent advances in RNA-based medicine have provided new opportunities for the global current challenge, i.e., the COVID-19 pandemic. Novel vaccines are based on a messenger RNA (mRNA) motif with a lipid nanoparticle (LNP) vector, consisting of high content of unique pH-sensitive ionizable lipids (ILs). Here we provide molecular insights into the role of the ILs and lipid mixtures used in current mRNA vaccines. We observed that the lipid mixtures adopted a nonlamellar organization, with ILs separating into a very disordered, pH-sensitive phase. We describe structural differences of the two ILs leading to their different congregation, with implications for the vaccine stability. Finally, as RNA interacts preferentially with IL-rich phases located at the regions with high curvature of lipid phase, local changes in RNA flexibility and base pairing are induced by lipids. A proper atomistic understanding of RNA-lipid interactions may enable rational tailoring of LNP composition for efficient RNA delivery.

Structure and drug binding of the SARS-CoV-2 envelope protein transmembrane domain in lipid bilayers.

An essential protein of the SARS-CoV-2 virus, the envelope protein E, forms a homopentameric cation channel that is important for virus pathogenicity. Here we report a 2.1-Å structure and the drug-binding site of E's transmembrane domain (ETM), determined using solid-state NMR spectroscopy. In lipid bilayers that mimic the endoplasmic reticulum-Golgi intermediate compartment (ERGIC) membrane, ETM forms a five-helix bundle surrounding a narrow pore. The protein deviates from the ideal α-helical geometry due to three phenylalanine residues, which stack within each helix and between helices. Together with valine and leucine interdigitation, these cause a dehydrated pore compared with the viroporins of influenza viruses and HIV. Hexamethylene amiloride binds the polar amino-terminal lumen, whereas acidic pH affects the carboxy-terminal conformation. Thus, the N- and C-terminal halves of this bipartite channel may interact with other viral and host proteins semi-independently. The structure sets the stage for designing E inhibitors as antiviral drugs.

Recombinant SARS-CoV-2 envelope protein traffics to the trans-Golgi network following amphipol-mediated delivery into human cells.

The severe acute respiratory syndrome coronavirus 2 envelope protein (S2-E) is a conserved membrane protein that is important for coronavirus (CoV) assembly and budding. Here, we describe the recombinant expression and purification of S2-E in amphipol-class amphipathic polymer solutions, which solubilize and stabilize membrane proteins, but do not disrupt membranes. We found that amphipol delivery of S2-E to preformed planar bilayers results in spontaneous membrane integration and formation of viroporin cation channels. Amphipol delivery of the S2-E protein to human cells results in plasma membrane integration, followed by retrograde trafficking to the trans-Golgi network and accumulation in swollen perinuclear lysosomal-associated membrane protein 1-positive vesicles, likely lysosomes. CoV envelope proteins have previously been proposed to manipulate the luminal pH of the trans-Golgi network, which serves as an accumulation station for progeny CoV particles prior to cellular egress via lysosomes. Delivery of S2-E to cells will enable chemical biological approaches for future studies of severe acute respiratory syndrome coronavirus 2 pathogenesis and possibly even development of "Trojan horse" antiviral therapies. Finally, this work also establishes a paradigm for amphipol-mediated delivery of membrane proteins to cells.

Microscopic Interactions of Melatonin, Serotonin and Tryptophan with Zwitterionic Phospholipid Membranes.

The interactions at the atomic level between small molecules and the main components of cellular plasma membranes are crucial for elucidating the mechanisms allowing for the entrance of such small species inside the cell. We have performed molecular dynamics and metadynamics simulations of tryptophan, serotonin, and melatonin at the interface of zwitterionic phospholipid bilayers. In this work, we will review recent computer simulation developments and report microscopic properties, such as the area per lipid and thickness of the membranes, atomic radial distribution functions, angular orientations, and free energy landscapes of small molecule binding to the membrane. Cholesterol affects the behaviour of the small molecules, which are mainly buried in the interfacial regions. We have observed a competition between the binding of small molecules to phospholipids and cholesterol through lipidic hydrogen-bonds. Free energy barriers that are associated to translational and orientational changes of melatonin have been found to be between 10-20 kJ/mol for distances of 1 nm between melatonin and the center of the membrane. Corresponding barriers for tryptophan and serotonin that are obtained from reversible work methods are of the order of 10 kJ/mol and reveal strong hydrogen bonding between such species and specific phospholipid sites. The diffusion of tryptophan and melatonin is of the order of 10-7 cm2/s for the cholesterol-free and cholesterol-rich setups.

On Drug-Membrane Permeability of Antivirals for SARS-CoV-2.

One of the key parameters required to identify effective drugs is membrane permeability, as a compound intended for an intracellular target with poor permeability will have low efficacy. In this paper, we leverage a computational approach recently developed by our group to study the interactions between nanoparticles and mammalian membranes to study the time of entry of a variety of drugs into the viral envelope of coronavirus as well as cellular organelles. Using a combination of all-atoms molecular dynamics simulations and statistical analysis, we consider both drug characteristics and membrane properties to determine the behavior of 79 drugs and their interactions with the viral envelope, composed of the membrane and spike protein, as well as five other membranes that correspond to various mammalian compartments (lysosome, plasma, Golgi, mitochondrial, and endoplasmic reticulum membranes). The results highlight important trends that can be exploited for drug design, from the relatively high permeability of the viral envelope and the effect of transmembrane proteins, to the differences in permeability between organelles. When compared with bioavailability data present in the literature, the model results suggest a negative correlation between time of permeation and bioavailability of promising drugs. The method is general and flexible and can be employed for a variety of molecules, from small drugs to small nanoparticles, as well to a variety of biological membranes. Overall, the results indicate that this model can contribute to the identification of successful drugs as it predicts the ability of compounds to reach both intended and unintended intracellular targets.

SARS-CoV-2 3a expression, purification, and reconstitution into lipid nanodiscs.

The SARS-CoV-2 3a protein is a putative ion channel implicated in virus life cycle and pathogenesis. We recently expressed, purified, and reconstituted 3a into lipid nanodiscs to solve its structure by cryo-EM to 2.1Å resolution. In this chapter, we describe methods we developed in order to facilitate the study of this protein in other laboratories. We emphasize factors that enabled rapid progression from gene sequence to reconstituted protein (3 weeks in the case of 3a) and provide general observations and tips for adapting these protocols to other membrane proteins of interest.

How Alcoholic Disinfectants Affect Coronavirus Model Membranes: Membrane Fluidity, Permeability, and Disintegration.

Atomistic molecular dynamics simulations have been carried out with a view to investigating the stability of the SARS-CoV-2 exterior membrane with respect to two common disinfectants, namely, aqueous solutions of ethanol and n-propanol. We used dipalmitoylphosphatidylcholine (DPPC) as a model membrane material and did simulations on both gel and liquid crystalline phases of membrane surrounded by aqueous solutions of varying alcohol concentrations (up to 17.5 mol %). While a moderate effect of alcohol on the gel phase of membrane is observed, its liquid crystalline phase is shown to be influenced dramatically by either alcohol. Our results show that aqueous solutions of only 5 and 10 mol % alcohol already have significant weakening effects on the membrane. The effects of n-propanol are always stronger than those of ethanol. The membrane changes its structure, when exposed to disinfectant solutions; uptake of alcohol causes it to swell laterally but to shrink vertically. At the same time, the orientational order of lipid tails decreases significantly. Metadynamics and grand-canonical ensemble simulations were done to calculate the free-energy profiles for permeation of alcohol and alcohol/water solubility in the DPPC. We found that the free-energy barrier to permeation of the DPPC liquid crystalline phase by all permeants is significantly lowered by alcohol uptake. At a disinfectant concentration of 10 mol %, it becomes insignificant enough to allow almost free passage of the disinfectant to the inside of the virus to cause damage there. It should be noted that the disinfectant also causes the barrier for water permeation to drop. Furthermore, the shrinking of the membrane thickness shortens the gap needed to be crossed by penetrants from outside the virus into its core. The lateral swelling also increases the average distance between head groups, which is a secondary barrier to membrane penetration, and hence further increases the penetration by disinfectants. At alcohol concentrations in the disinfectant solution above 15 mol %, we reliably observe disintegration of the DPPC membrane in its liquid crystalline phase.